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ABSTRACT 

Accurate galaxy stellar masses are crucial to better understand the physical mecha- 
nisms driving the galaxy formation process. We use synthetic star formation and metal 
enrichment histories predicted by the GALFORM galaxy formation model to investigate 
the precision with which various colours (m a — nib) can alone be used as diagnostics 
of the stellar mass-to-light ratio. As an example, we find that, at z = 0, the intrinsic 
(Bj435 w — V/606tu) colour can be used to determine the intrinsic rest-frame T^-band 
stellar mass-to-light ratio (log 10 IV = l°gio[(M /Mq)/(Lv/Lvq)]) wrt h a precision of 
uigY — 0.06 when the initial mass function and redshift are known beforehand. While 
the presence of dust, assuming a universal attenuation curve, can have a systematic 
effect on the inferred mass-to-light ratio using a single-colour relation, this is typically 
small as it is often possible to choose a colour for which the dust reddening vector 
is approximately aligned with the (m a — nib) — log 10 IV relation. The precision with 
which the stellar mass-to-light ratio can be recovered using a single colour diagnostic 
rivals implementations of SED fitting using more information but in which simple pa- 
rameterisations of the star formation and metal enrichment histories are assumed. To 
facilitate the wide use of these relations, we provide the optimal observer frame colour 
to estimate the stellar mass-to-light ratio, along with the associated parameters, as a 
function of redshift (0 < z < 1.5) for two sets of commonly used filters. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high- redshift - ultraviolet: galaxies 
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1 INTRODUCTION 

The accurate determination of the stellar masses of galaxies 
at different redshifts is crucial to improve our understanding 
of the process of galaxy formation and evolution. 

The physical properties of galaxies, including the stellar 
mass-to-light ratio, are most often determined using spectral 
energy distribution (SED) fitting (see Walcher et al. 2011 
for a recent overview), using broad-band photometry (e.g. 
Taylor et al. 2011), spectroscopy (e.g. Tojeiro et al. 2007, 
Pacifici et al. 2012), or spectral indicies (e.g. Kauffmann et 
al. 2003). The determination of the various properties (in- 
cluding the stellar mass-to-light ratio) involves the compar- 
ison of the observed SED with a library of templates. This 
comparison has been done either using a single best-fitting 
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template (e.g. the template/model which minimises \ 2 ) or 
carrying forward the full posterior probability distribution 
which contains all the information inferred about the stellar 
mass-to-light ratio (e.g. Kauffmann et al. 2003, Taylor et al. 
2011.) 

Each template/model is constructed using stellar popu- 
lation synthesis (SPS) modelling (e.g. Tinsley & Gunn 1976, 
Bruzual & Chariot 1993). Stellar evolution models are com- 
bined with a choice of initial mass function (IMF) to produce 
synthetic single-age stellar population (SSP) SED models as 
a function of age and metallicity (e.g. Bruzual & Chariot 
2003, Fioc &: Rocca-Volmerange 1997). By combining the 
SSP models with star formation and metal enrichment his- 
tories the intrinsic SED of a composite stellar population 
(CSP) can be assembled. By including absorption and emis- 
sion due to dust and gas in the ISM and IGM (e.g. Madau 
et al. 1995) realistic SEDs of galaxies can be constructed. In 
its most basic incarnation, SED fitting assumes a template 
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library built using a simple parameterisation of both the star 
formation and metal enrichment histories with some param- 
eter^) describing the effect of dust attenuation. The choice 
of star formation history (SFH) and metal enrichment his- 
tory parameterisations varies between studies depending on 
the available photometry and/or spectroscopy. For example, 
at very-high redshift (z > 5) where there is little high signal- 
to-noise photometry, a constant SFH is sometimes assumed 
(e.g. Gonzalez et al. 2011). However, a more common pa- 
rameterisation is a decaying exponential: oc e _t ' T (which has 
two free-parameters, the age t and characteristic timescale 
r) which is able to account for galaxies with little residual 
star formation. Further improvements include the addition 
of bursts (e.g. Kauffmann et al. 2003, Brinchmann et al. 
2004, Salim et al. 2007, da Cunha et al. 2008, Walcher et 
al. 2009) or the truncation of star formation entirely. Fur- 
ther, some studies (e.g. Tojeiro et al. 2007) avoid the need 
for SFH and metal enrichment parameterisations altogether 
and instead construct CSPs from a series of discrete star 
formation episodes. 

A different alternative to assuming a parameterisation 
is to build a training set (i.e. a set of observations with 
known properties). While it is possible to build an empirical 
training set in some cases (for example when the aim is 
to determine the redshift) in practice this is not possible 
for the stellar mass-to-light ratio as this cannot be directly 
measured (unlike the redshift). However, we can instead use 
galaxy formation models to produce physically motivated 
star formation and metal enrichment histories that can be 
used to generate what is essentially an artificial training set. 
This approach is of course only effective if the simulation is 
a good description of the real Universe. 

While it is possible to use an artificial training set as 
the basis of full SED fitting (i.e. using all available informa- 
tion) in this paper we focus on deriving a simple relationship 
between a single-colour and the stellar mass-to-light ratio. 
This is particularly useful for spectroscopic surveys of galax- 
ies which may not have extensive multiband photometry. 
In addition, mass estimates based on a single-colour rela- 
tion (unlike SED fitting with more than 2-bands) are easily 
implemented by other studies permitting a straightforward 
comparison between samples. However, any empirical or the- 
oretically defined single-colour - mass-to-light ratio relation 
is based on a specific population and may not be suitable 
for populations defined using an alternative criteria. 

The use of a single-colour diagnostic has previously been 
considered in literature (at low redshift) and the strength of 
this approach has been demonstrated by, for example, Bell 
& de Jong (2001) and more recently Taylor et al. (2011). 
For example, Taylor et al. (2011) compared stellar mass- 
to-light ratios derived from SED fitting using broad-band 
optical photometry collected by the Galaxy and Mass As- 
sembly (GAMA, Driver et al. 2009, 2011) survey with the 
rest-frame (g — i) colours finding a strong correlation. 

It is also important to point out that stellar masses 
derived from any photometric method are sensitive to the 
choice of initial mass function (IMF). For example, stellar 
masses determined assuming a Salpeter (1955) IMF are, due 
to enhanced contribution of faint low-mass < 0.5 Mq stars, 
typically around 0.2 dex greater than those derived assuming 
a Chabrier (2003) IMF. This problem is exacerbated further 
if the IMF is not universal as suggested by several different 



studies (e.g. Hoversten & Glazebrook 2008, van Dokkum 
2008, Wilkins et al. 2008a, Lee et al. 2009, Cappellari et al. 
2012). Stellar mass estimates are also sensitive to uncertain- 
ties in our understanding of stellar evolution (e.g. Conroy et 
al. 2009) and thus vary depending on the choice of stellar 
population synthesis code (e.g. Maraston et al. 2003, Tonini 
et al. 2010). 

In this work we are interested firstly in determining the 
accuracy and precision with which a single colour can be 
used to infer the stellar mass-to-light ratio, the theory of 
which is described in Section [2] To achieve this we use the 
GALFORM galaxy formation model to generate realistic star 
formation and metal enrichment histories by following the 
physical processes that affect the baryonic component of the 
Universe. These are in turn used to produce synthetic intrin- 
sic photometry for a large number of galaxies. The GALFORM 
model and filter sets utilised are described in §3.1| and §3.2| 
respectively. In Section HI we demonstrate the precision with 
which a carefully chosen intrinsic colour can be used to de- 
termine the stellar mass-to-light ratio. We then investigate 
the effect of dust in §4.2| and in Section [5] we determine, and 
present (along with the associated parameters), the optimal 
choice of colour with which to probe the mass-to-light ratio 
over z = 0.0 — > 1.5. In Section [6] we compare the expected 
precision of a single-colour based diagnostic of the mass-to- 
light ratio with recent estimates of the accuracy and uncer- 
tainties of SED fitting. Finally in [7] we present our conclu- 
sions. Throughout this work we employ the AB magnitude 
system (Oke & Gunn 1983) and assume a Salpeter (1955) 
IMF, i.e. £(m) = AN /dm oc m" 2 ' 35 . 



2 THE REST-FRAME OPTICAL COLOUR AS 
A DIAGNOSTIC OF THE MASS-TO-LIGHT 
RATIO 

The mass-to-light ratio of a stellar population is driven pre- 
dominantly by the distribution of stellar masses and ages; 
the specific (i.e. per unit mass) luminosity of high-mass stars 
is much larger than for low-mass stars. The distribution of 
stellar masses is in turn driven by the star formation his- 
tory (SFH) and initial mass function. Galaxies with pro- 
tracted but declining star formation histories will contain, 
relative to galaxies with substantial recent star formation, 
more low-mass, lower-luminosity stars due to their longer 
main-sequence lifetimes. 

Fig. n] shows the specific SED (per unit stellar mass) 
assuming an exponentially decaying star formation history 
(oc e" t/T ) with r = lGyr and log 10 (i/Gyr) = -1.0 -> 1.0 
in increments of Alog 10 (t/Gyr) = 0.2. As the stellar pop- 
ulation ages the specific SED decreases (the stellar mass- 
to-light ratio increases) and the rest-frame optical colour 
reddens. The correlation between stellar mass-to-light ratio 
and rest-frame optical colour can be seen more clearly in 
Fig. [5] where evolutionary tracks (f/Gyr = — > t unl (z = 
0) = 13.5) in the rest-frame HST ACS {B fi35w - U/eoe™) 
colour - mass-to-light ratio plane are shown for various ex- 
ponentially decaying star formation histories {ip oc e~*' T 
with r/Gyr € {0.1,0.3, 1.0,3.0, 10.0}) and two metallicities 
{Z/Zq = 1 and Z/Zq — 0.02) constructed using the pe- 
GASE.2 SPS model (Fioc & Rocca-Volmerange 1997, 1999) 
and assuming a Salpeter IMF. 
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Figure 1. Upper panel: The specific SED (i.e. per unit current stellar mass) assuming an exponentially decaying star formation history 
(<x e~*/ T ) with t = 1 Gyr and log 10 (t/Gyr) = —1.0 — > 1.0 in increments of A log 10 (t/Gyr) = 0.2 (plotted using black to light grey lines). 
The right-hand axis shows the corresponding mass-to-light ratio (normalised by the V-band solar luminosity and solar mass). Lower 
panel: Transmission functions of the HST ACS and WFC3 filters considered in this study. 



For stellar populations with colours (.8/435™ — V/606™)> 
(rjthe colour and mass-to-light ratio are correlated, i.e. for 
a given colour the range of stellar mass-to-light ratios is 
substantially smaller than the full range of mass-to-light ra- 
tios. While stellar populations with lower metallicities ex- 
tend over a smaller range of mass-to-light ratios and colours, 
the mass-to-light ratio at the same colour is roughly similar, 
suggesting that using an optical colour to infer the mass-to- 
light is somewhat insensitive to the metal enrichment his- 
tory. 

However, it is difficult to derive an accurate single- 
colour mass-to-light ratio relation, together with an estimate 
of the uncertainty, from such simple and unrealistic star for- 
mation and metal enrichment histories. Instead it is more 
useful to use a galaxy formation model to produce a realis- 
tic library of galaxy star formation histories. 



3 MODELLING APPROACH 

3.1 Galaxy Formation Model 

We predict the stellar masses and colours of galaxies at dif- 
ferent redshifts in a ACDM universe using the GALFORM 
semi-analytical galaxy formation model initially developed 

1 Colours < typically occur only for very young stellar popula- 
tions. 
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Figure 2. Evolutionary tracks (t = — > t un [(z = 0)) of stel- 
lar populations in the rest-frame (/435ui — /606«i) - V-band 
mass-to-light ratio plane, created assuming a Salpeter (1955) 
IMF, utilising the PEGASE.2 SPS model and assuming several ex- 
ponentially declining star formation histories (ip ex e~ t ' T with 
-r/Gyr e {0.1, 0.3, 1.0, 3.0, 10.0}) and two metallicities (Z/Z Q = 1 
and Z/Zq = 0.02) as indicated by the key. 



by Cole et al. (2000). Semi-analytical models solve physi- 
cally motivated differential equations to follow the fate of 
baryons in a universe in which structure grows hierarchi- 
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cally through gravitational instability (see Baugh 2006 for 
an overview of hierarchical galaxy formation models). 

GALFORM follows the main processes which shape the 
formation and evolution of galaxies. These include: (i) the 
collapse and merging of dark matter haloes; (ii) the shock- 
heating and radiative cooling of gas inside dark matter 
haloes, leading to the formation of galaxy discs; (iii) qui- 
escent star formation in galaxy discs; (iv) feedback from su- 
pernovae, from active galactic nuclei (AGN) and from pho- 
toionization of the intergalactic medium (IGM); (v) chem- 
ical enrichment of the stars and gas; (vi) galaxy mergers 
driven by dynamical friction within common dark matter 
haloes, leading to the formation of stellar spheroids, which 
also may trigger bursts of star formation; (vii) bursts of star 
formation in dynamically unstable discs. The end product 
of the calculations is a prediction for the number and prop- 
erties of galaxies that reside within dark matter haloes of 
different masses. These properties include the star forma- 
tion and metal enrichment histories. 

In this paper we focus our attention on the Bower et al. 
(2006, B06) variant of GALFORM. Key features of the B06 
model include (i) a time scale for quiescent star formation 
that scales with the dynamical time of the disk and which 
therefore changes significantly with redshift (see Lagos et 
al. 2011 for a study of different star formation laws in qui- 
escent galaxies), (ii) bursts of star formation occur due to 
both galaxy mergers and when disks become dynamically 
unstable, and (iii) negative feedback from both supernovae 
and AGN (see Benson et al. 2003 for a discussion of the 
impact of feedback has on the luminosity function of galax- 
ies) . The onset of the AGN suppression of the cooling flow 
is governed by a comparison of the cooling time of the gas 
with the free-fall time for the gas to reach the centre of 
the halo. Cooling is suppressed in quasi-static hot haloes 
if the luminosity released by material accreted onto a cen- 
tral supermassive black hole balances the cooling luminosity 
(see Fanidakis et al. 2011 for a full description of black hole 
growth in the model). Bower et al. (2006) adopt the cosmo- 
logical parameters of the Millennium Simulation (Springel 
et al. 2005), which are in broad agreement with constraints 
from measurements of the cosmic microwave background ra- 
diation and large scale galaxy clustering (e.g. Sanchez et al. 
2009): fi = 0.25, A = 0.75, Q b = 0.045, erg = 0.9 and 
h = 0.73, such that the Hubble constant today is Ho = 100 h 
kms _1 Mpc _1 . The B06 model parameters were fixed with 
reference to a subset of the available observations of galaxies, 
mostly at low redshift (see Bower et al. 2010 for a discussion 
of parameter fitting). This model successfully reproduces the 
galaxy stellar mass function up to z — 4.5 and the proper- 
ties of red galaxies at 0.7 < z < 2.5 (Gonzalez-Perez et al. 
2009, 2011). For further details we refer the reader to B06. 

For the work presented here, we have implemented the 
B06 model using merger trees generated with the Monte 
Carlo algorithm introduced by Parkinson et al. (2008). We 
have also changed the default SSP used to produce the SEDs 
of galaxies, using here PEGASE.2 combined with a Salpeter 
(1955) IMrQ A Salpeter IMF is chosen over the default 
(Kennicutt 1983) IMF due to its wide use in the literature. 



These changes have only a minor effect on the intrinsic dis- 
tribution of colours and do not affect the optimum choice 
of colour with which to probe the stellar mass-to-light ra- 
tio. Assuming a Salpeter IMF however increases the stellar 
mass-to-light ratios (relative to a Kennicutt IMF) and thus 
affects the normalisation of any relation between colour and 
stellar mass-to-light ratio. The adoption of Salpeter IMF 
modifies the predicted luminosity function in the expected 
way, introducing a shift of about 0.4 mag in the 6-band lu- 
minosity function at z=0 with respect to observations. 

In this work, we uniquely use the unextincted magni- 
tudes output by GALFORM. The attenuation of starlight by 
dust is modelled in GALFORM in a physically self-consistent 
way, based on the results of a radiative transfer calculation 
for a realistic geometry and predicted galaxy sizes and dust 
masses (Cole et al. 2000, Lacey et al. 2011, Gonzalez-Perez et 
al. 2012). Nevertheless, here we do not make use of this part 
of the modelling (see Mitchell et al., in prep., for a discussion 
on the stellar mass estimate using the standard dust treat- 
ment in GALFORM). Instead, here we apply a stand-alone 
model for dust extinction using the Calzetti et al. (2000) 
reddening law (see [ 4.2 1, one of the most widely used tech- 
niques for analysing the dust content of observed galaxies. 
This allows us to demonstrate the effect that a universal 
attenuation curve has on the stellar mass estimation. 



3.1.1 Sample definition 

Throughout this work we consider a stellar mass limited 
sample, choosing \og x0 {M*/Mo) > 9/j This limit is chosen 
to roughly reflect the mass range accessible to typical HST 
surveys at z = 0.5 — 1.5. 

The merger trees used in the model are generated using 
a Monte Carlo approach, for which weights are given to the 
haloes quantifying their abundance based on the halo mass 
function. In turn, each galaxy has the associated weight of 
its host halo encoding the relative probability with which it 
would occur. Thus, a histogram of the stellar mass of each 
galaxy in the sample incorporating these weights will then 
reproduce the galaxy stellar mass function. We carry these 
weights through our entire analysis to ensure that the sample 
truly reflects a mass-limited criteria. 



3.1.2 The predicted distribution of mass-to-light ratios 

The distribution of intrinsic rest-frame V-band mass-to-light 
ratios (log 10 Fv) predicted by GALFORM is shown in Fig. [3] 
The distribution spans a range of A log 10 Fy — 0.61 between 
the 10 th to 90 th percentiles and has a standard deviation of 
°"iog 10 r = 0.26, though the distribution is non-gaussian and 
asymmetric. The bi-modality simply reflects the population 
of active and passive galaxies predicted by the model, cor- 
responding to 'red' and 'blue' galaxies. 



3.2 Filter Sets 

In this work we consider two sets of filters. The HST fil- 
ter set employs a combination of HST ACS and WFC3 



2 Note that B06 used the SSPs from Bruzual & Chariot (1999) 
combined with a Kennicutt (1983) IMF. 



3 In j j4.l| we also consider the effect on the recovered parameters 
of changing this mass cut finding only a small effect. 
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Figure 3. The distribution of intrinsic rest-frame V-band mass- 
to-light ratios (IV) for galaxies with log 10 (M / M®) > 9 at z = 
predicted by the Bower et al. (2006) version of the GALFORM 
galaxy formation model assuming a Salpeter IMF and the PE- 
GASE.2 SPS model. The vertical lines denote the 15.87 th , 50 th , 
and 84.13 th percentiles of the distribution from left to right re- 
spectively. 



filters (WFC3.UVIS: UV /225 », UV /275m , U /336 », ACS: 

B/435m, V/606m, if775w, ZfS50lp, WFC3.NIR: Y/i05u>, J/125tu, 

H/ieouj)- The choice of these filters is motivated by the exis- 
tence of publically availability deep imaging in all or a sig- 
nificant combination of these filters. The GAM A filter set 
is based on the observations accessible to the GAMA sur- 
vey and includes the Galaxy Evolution Explorer (GALEX, 
Martin et al. 2005) FUV and NUV, Sloan Digital Sky Survey 
(SDSS, Fukugita et al. 1996) ugriz and the United Kingdom 
Infrared Telescope ( UKIRT) Wide Field Camera (WFCAM, 
Casali et al. 2007) YJHK filters (Hewett et al. 2006). The 
filter transmission functions are shown in Fig. [T] 



4 SINGLE-COLOUR MASS-TO-LIGHT RATIO 
RELATION 

As an example we consider the usefulness of the intrinsic 
observed-frame (B/4351D — V/606w) colour as a diagnostic of 
the rest-frame V-band mass-to-light ratio at z = 0. The 
(B/435u,— V/606jd) colour is chosen to demonstrate the use- 
fulness of the technique as it represents the optimum choice 
of colour at z = (see Section pt with which to measure 
the mass-to-light ratio using the HST filter set. Comparing 
the intrinsic observed- frame (B/435™— V/606™) colour and 
the rest-frame V-band mass-to-light ratio (log 10 Fv) pre- 
dicted by GALFORM at Z = 0, as shown in Fig. HI reveals 
the two quantities are highly correlated (r ~ 0.9, where 
r is Pearson's sample correlation coefficien1r| . Using lin- 
ear regression, taking into account the abundance of each 
galaxy, we can fit the relationship between an arbitrary 
colour (m a — nib) and the V-band mass-to-light ratio by 
a linear relation, i.e. 



logio r i 



log w [(M/M Q )/(L v /L VQ )} 



(1) 



4 A value \r X y\ = 1.0 implies the relationship between the vari- 
ables x and y can be perfectly described by a linear equation, 
while r = 0.0 on the other hand implies there is no linear corre- 
lation between the variables. 



Figure 4. The relationship between the (1^435^ — Vj606™) colour 
and the V-band mass-to-light ratio (log 10 Ty) for synthetic galax- 
ies, with log 10 (M/Mq) > 9.0, extracted from GALFORM at z = 0. 
The vertical lines denote the 15.87 th -84.13 th percentile range of 
distribution in several colour bins. The solid (green) line is the 
best-fitting linear relation between the colour and mass-to-light 
ratio. The coloured curves show the predicted tracks for expo- 
nentially decreasing star formation histories as in Fig. [2] The 
greyscale points show the location of a representative sub-sample 
of galaxies. 



pi X (m a - m b ) + p 2 . 



(2) 



For the (-B/435U, — Vfeoew) colour / V-band mass-to-light ra- 
tio the parameters pi and p 2 are 1.10 and —0.50 respectively 
and the linear fit is displayed in Fig. [4] Reassuringly this lin- 
ear fit is coincident with the predictions made in Section [2] 
assuming simple exponentially decaying star formation his- 
tories (as shown in Fig. [2] and Fig. Hi. 

If we now consider the residual distribution (R, shown 
in Fig. pi between the mass-to-light ratio estimated from 
this linear relation and the true mass-to-light ratio, defined 
generally as, 



R = log 10 IV - (pi(m a - m b ) +P2), 



(3) 



we find that, for the (B/435™— V/606uj) colour, the standard 
deviation of R is an ~ 0.059. This is a factor ~ 4x smaller 
than the standard deviation of the stellar mass-to-light ratio 
distribution (o"i 9 r = 0.261, see j |3.1.2[ |. The residual distri- 
bution (as shown in Fig. pMI has an approximately gaussian 
shape (Fig.[5]also shows a normal distribution with the same 
mean and standard deviation as R) except for the tails of the 
distribution. This analysis suggests that a perfectly measured 
intrinsic colour can be used to estimate the mass-to-light ra- 
tio Yv with a precision of ~ 14%. An additional example is 
shown in Fig. [6] where a similar relation is evident at z = 1. 



4.1 Sensitivity to stellar mass cut 



As noted in §3.1.1| our sample is effectively mass limited, 
with log lo (M*/M0) > 9. By changing this criteria we can 
assess the sensitivity of the parameters p\ and p 2 - We find 
that both pi and p 2 vary by only a small amount (pi = 
1.10 -» 1.07, p 2 = -0.50 -¥ -0.52) as the mass limit is 
increased from log lo (M*/M0) > 9 — > 11. 
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Figure 5. The distribution of residuals (Ri), that is, the differ- 
ence between the true mass-to-light ratio T and that inferred from 
the intrinsic colour using the best-fit linear relation, Eq. [I] The 
solid grey curve represents a normal distribution with the same 
mean, standard deviation and area as Ri. The left and right ver- 
tical dashed lines denote the 15.87 th -84.13 th percentiles respec- 
tively, while the dotted line denotes the median. 
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Figure 6. Similar to Fig.|4]but showing instead the relationship 
between the (z/ssoip — Jfi25w) ar >d the V-band mass-to-light ra- 
tio (log 1( )rV) for synthetic galaxies, with log lo (M/M0) > 9.0, 
extracted from GALFORM at z = 1. 



4.2 Dust 

In the preceding analysis we demonstrated that the intrinsic 
colour of a galaxy provides a robust diagnostic of the stellar 
mass-to-light ratio. However, the observed SEDs of galaxies, 
and thus their colours and observed mass-to-light ratios are 
altered by dust, potentially affecting the reliability of a single 
colour as a diagnostic. Here we analyze the effects of dust 
assuming an attenuation law with a universal shape, with 
the overall amount of dust extinction specified by a single 
parameter which varies from galaxy to galaxy. 

Over the UV, optical and NIR, dust acts to preferen- 
tially absorb light at shorter wavelengths. One effect of dust 
is then to redden the observed SED and thus the colour rel- 
ative to the intrinsic, i.e. (m a — m b ) ohB = {m a — rnb) lnt + 
(A a — Ab), where (A a — At) > 0.0. However dust, by ab- 
sorbing some fraction of the light, also acts to increase the 
observed mass-to-light ratio relative to the intrinsic value, 



(i.e. r° bs > T int ). Thus the effect of dust is to shift an indi- 
vidual galaxy in the colour-F ratio plane by a vector d with 
orientation < 9 < n/2 (with respect to the log 10 F axis). 
That is, dust acts in approximately the same direction as the 
(m a — mf,)-log 10 Tv relation. 

The direction of this vector can be estimated analyti- 
cally as follows: the relationship between the observed and 
intrinsic mass-to-light ratio is, 

log 10 (r£ bs ) = Iog 10 (Ti? t )+0.4£ fl _yfcv, (4) 



where Eb-v is the colour excess and ki is the effective red- 
dening for a filter i such that ki = Ai/ Eb-v- The extinction 
Ai is defined as 



a obs int 

Ai — rrii — rrii 

which is in turn given by, 

obs ini 

rrii — rrii 
where, 



-2.5 log 10 Q, 



Q 



i: , n -0AE B _ v k x df 



T/™10 



Tf: 



V 



(5) 



(6) 



(7) 



where T is the filter response curve and fl nt is the intrinsic 
spectrum of the sourcerl 

The relationship between between the observed and in- 
trinsic colour is, 



(m a — m b y 



Thus, the dust vector is then given by, 
d = (E B -v(ka-kb),0AE B -vkv) 



{m a 


- m 6 ) mt + {Aa - A b ) 


(8) 


(m a 


- m b ) lnt + E B -v(k a - fcfc). 


(9) 



which has a gradient, 
di = 0.4 x k v /{k a - k b ). 



(10) 



(11) 



Using this formalism, and assuming a Calzetti et al. 
(2000) attenuation curvaj we determine that the gradient 
of the dust vector d, for the (B/435™— V/606™) colour at 
z = 0, is di = 1.30 (c.f. pi = 1.10), i.e. while the (m a — m b )- 
log 10 Tv relation and dust vector point in approximately the 
same direction they are not perfectly aligned. Qualitatively 
the presence of a dust vector not perfectly aligned with the 
relation will then have two effects: it will both increase the 
scatter in the residual distribution (i.e. decrease the preci- 
sion of the {m a — mj,)-log 10 Fv relation) and introduce a 
systematic shift. 

The direction of the systematic shift depends on the 



5 Examination of Eq.|7|shows that ki, for a given filter, is sensitive 
to the choice of Eb—v & n d /i nt and therefore differs for every 
distinct /" lt; . However, over the range Eb—v = 0-0 — > 2.0 and for 
fv oc u^ +2 , with ft = —2 — > 2, ki varies only at the percent level. 
A notable exception where this assumption breaks down is when 
the light in a given filter is dominated by emission from a small 
wavelength interval; for example when the equivalent width of an 
emission line is comparable to the width of the filter, or where 
the filter itself stradles a spectral break. 

6 The shapes of the commonly used dust attenuation curves, in- 
cluding the Calzetti et al. (2000) starburst, Milky Way, SMC and 
LMC (e.g. Pei et al. 1992) curves, are all similar in the rest-frame 
optical where the optimal colour is typically found. 
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angle between the dust vector gradient di and pi. For ex- 
ample, if the gradient of the dust vector is smaller than that 
of the intrinsic relation (i.e. di < pi) the mass-to-light ratios 
inferred using the intinsic relation will be systematically too 
high, i.e. the stellar masses will be overestimated. The size of 
the systematic shift (51og 10 IV) can be derived analytically, 



<51og 10 r> = (k a - h)(di -pi)E B - 



d 2 x E B - 



(12) 



For the (B/4 35u ,— V/606m) relation this gives d 2 = 0.26, sug- 
gesting only a small (< 0.1) systematic shift for typical val- 
ues of Eb-v (< 0.3). However, for extremely dusty galaxies 
the mass-to-light ratio will be overestimated. It is worth not- 
ing that the misalignment between d and the single-colour-r 
relation in this case is in fact greater than every other opti- 
mum choice of colour (see the proceeding section) . 

4-2.1 The effect of normally distributed dust 

The effect of dust can be further illustrated by applying dust 
extinction to each galaxy according to some distribution of 
Eb-v- For this we assume each value of Eb-v is drawn from 
a truncated normal distribution such that Eb-v > 0.0, i.e. 
Vv(Eb-v) =N(p,cr 2 )[E B -v > 0Q Fig. [7] demonstrates the 
effect of assuming this form of dust distribution on the resid- 
ual distribution (for the T5f43s w — ~Vf606w relation at z = 0) 
assuming various values of /i (G {0.0,0.1,0.25, 0.5}tj and a 
(e {0.1,0.25}). Increasing 11 = — > 0.5 causes a systematic 
increase in R of 0.13 (the same as the value predicted using 
Eq. 12 1 while increasing the scatter a = — > 0.25 results 
in the standard deviation of R increasing from w 0.06 to 
r;0.09. 



4.3 Photometric uncertainties 

Another advantage of a linear relation between the colour 
(m a — mb) and the stellar mass-to-light ratio log 10 Tv is that 
it is easy to show that log 10 Tv is gaussian assuming both 
the distribution of R and the noise on the colour are also 
gaussian. The mean of loginJV is clearly pi(m a — m,b) + p 2 
and the uncertainty is then 



2 
o-i a r 



22.2 

Vl°~ab + PR, 



(13) 



where o a b is the uncertainty on the colour (m a — mb) 



5 CHOOSING THE OPTIMAL COLOUR 

We have now shown that even in the presence of dust a 
single observed colour provides a robust diagnostic of the 
stellar mass-to-light ratio. It is useful then to determine, for 
a series of redshifts (0 < z < 1.5 for the HST filter set and 
< z < 0.5 for the GAMA filter set), the optimum colour 
with which to estimate the mass-to-light ratio together with 
the associated parameters. When determining the optimum 



7 Here Af(fJ., cr 2 ) denotes a normal distribution with mean fi and 
variance <r 2 while the term [Eb—v > 0] forces the probability on 
un- physical negative values of Eb — v to be zero using the Iverson 
bracket notation. 

8 Because we have truncated the distribution fi is no longer equal 
to the mean/ median value of Eb-v ■ 

9 This can be shown by integrating over the colour distribution. 



Pr(E(B-V))=N(0.5,0.25 2 )[E(B-V)>0] 
meon(R)=0.1 33 
std(R) =0.086 
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Pr(E(B-V))=N(0.25.0.1 2 )[E(B-V)>0] 
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sld(R) =0.064 
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mean(R)=0.03 
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meon(R)=0.049 

std(R)=0.07 
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J u_ 




-0.4 -0.2 0.2 

R =| og,o(r v , obs )-(p,+p 2 x(B f4 3 5 „-' 



0.4 



Figure 7. The resulting distribution of stellar mass-to-light ratio 
residuals R (Eq.pl after applying dust using Eb—v drawn from a 
truncated normal distribution (Pt(Eb-v) = A/"(/i, <j 2 )[Eb—v ^ 
0]) for several different values of /j, (e {0.0,0.1,0.25,0.5}) and 
o" (6 {0.1,0.25}). In all panels the shaded histogram shows the 
intrinsic distribution of residuals (as also shown in Fig. pi. 



colour we consider all 

in each of the filter sets described in 93.2 



the possible combinations of filters 
(i.e. 



we consider 

45 possible colours based on the 9-filter HST set and 55 
colours from the 10-filter GAMA set). Colours probing the 
rest-frame UV typically provide the worst correspondence 
with the mass-to-light ratio (R > 0.13, and in the case of 
the FUV ~ NUV colour at z = R ~ 0.3). This is predom- 
inantly because the UV is strongly affected by the recent 
star formation history of the galaxy which does not neces- 
sarily correlate with overall mass-to-light ratio. At z > 0.3 
the FUV (and at z > 1.1 the UV/225tu) band is also affected 
by the Lyman-limit break. 

In all redshift cases a range of filter combinations pro- 
duce roughly similar residual distributions, with or < 0.07. 
For example, at z = both the (B/435 m — V/606u>) and 
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(B/435 m — if775w) yield or < 0.06. Because of this we can 
also dictate a preference for colours for which the difference 
between the (m a — mb)-log 10 Fv relation slope and the gra- 
dient of the dust vector (di) is minimised, thereby reducing 
the systematic shifts expected due to dust. 

For both filter sets and a range of redshifts (GAM A: 
Z = 0.0- 0.5, HST: z = 0.0 - 1.5) Table [l] gives the best- 
fitting colour and gives the parameters pi, p%, d\ and cfe 
for each redshift. In virtually all cases the optimal choice 
of colour covers the rest-frame 0.4 — 0.7/im range. In each 
case the standard deviation of the relevant intrinsic residual 
distribution is or < 0.07. In most cases the dust vector is 
closely aligned (typically \d%\ < 0.15) with the (m a — mt)- 
log 10 IV relation with the greatest (fractional) offset being 
the (B/435„— V/606m) colour at z — (considered in the 
previous section). 

Thus far we have derived the stellar mass-to-light ratio 
- single colour relation using the rest-frame V-band lumi- 
nosity. To determine the stellar mass it is then necessary 
to apply a fc-correction to an observed luminosity (to de- 
termine the rest-frame 1^-band luminosity) introducing an 
additional source of uncertainty. However, we can instead 
define the mass-to-light ratio - single colour relation using 
an observed frame luminosity thereby removing the need for 
a fc-correction. This is detailed in Appendix [A"| 



6 COMPARISON WITH OTHER STUDIES 

6.1 Other determinations of single-colour 
relations 

It is interesting to assess whether our relation(s) are simi- 
lar to those derived empirically in the literature. Taylor et 
al. (2011, hereafter Til) define an empirical relation using 
stellar masses and (fc-corrected) colours from the Galaxy 
and Mass Assembly (GAMA) spectroscopic survey (Driver 
et al. 2011). The GAMA survey currently incorporates UV, 
optical and near-IR photometry from GALEX, SDSS and 
UKIDSS respectively with spectroscopy coming from a va- 
riety of sources (predominantly from SDSS and AAO). Til 
derives stellar masses using a Bayesian SED fitting method- 
ology using only optical (SDSS ugriz) photometry. The Til 
single-colour - mass-to-light ratio relation is defined using 
the rest-frame (g — i) colour and i-band luminosity such 
that, 



log 10 V^ ab = 0.7 x (g - i) r + 1.15, 



(14) 



where log 10 F^f, is defined as the mass-to-light ratio using 
AB centric units (i.e. luminosity has units of Li(Mt — 0.0)). 
At 2 = we find that the (g — r) colour is favoured over 
(g — i) due to the slightly smaller value of R. However, we 
can nevertheless easily define a relation between the (g — i) 
colour and and i-band luminosity (at 2 = 0), finding: 



log 10 r,, ab = 0.80 X (g - i) r + 1.41, 



(15) 



where log 10 Ti ia b is defined using the same AB centric unit 
system in Eq. |14| 

Clearly there is some disagreement between this relation 
and the Til relation. However there are several important 
differences in their construction which need to be addressed. 
Firstly, while our relation is defined using a sample of galax- 
ies at 2 = (such that the observed g — i colour is also the 



rest-frame g — i colour) the Til (GAMA) sample encom- 
passes a range of redshifts (z < 0.65; median z = 0.2). If 
instead we consider the rest-frame (g — i) colour at the me- 
dian GAMA redshift (z = 0.2) we find, 



logio r i,*b = °- 74 x (9 - i)r + 1-41, 



(16) 



the slope of which is now consistent with that derived by 
Til. Secondly, stellar masses in Til were derived assuming 
a Chabrier (2003) initial mass function (IMF) while, in our 
case, GALFORM was run assuming a Salpeter (1955) IMF. As 
noted in the introduction stellar masses estimated assuming 
a Salpeter IMF are typically around 0.2 dex larger than those 
assuming the Chabrier (2003) IMF. Making this correction 
to the normalisation of the Til relation (i.e. 1.15 + 0.2 = 
1.35) then yields good agreement between this single colour 
mass-to- light ratio relation and that of Til. 



6.2 Multi-wavelength broadband SED fitting 

It is also useful to comment on the precision of a single- 
colour diagnostic compared with typical broad-band SED 
fitting techniques. Several recent studies have attempted 
to determine the accuracy and precision (uncertainty) with 
which broad-band SED modelling can recover various phys- 
ical properties, including the stellar mass-to-light ratio (e.g. 
Pacifici et al. 2012, Pforr et al. 2012, Wilkins et al. in prep). 

Wilkins et al. in prep uses the same library of syn- 
thetic galaxies considered in this paper to investigate the 
the accuracy and precision (uncertainty) with which the 
stellar mass-to-light ratio can be recovered for various ob- 
servational factors (including the availability of a priori red- 
shifts, rest-frame SED coverage and S/N), modelling choices 
(including star formation history parameterisation, statisti- 
cal indicator) and as a function of various galaxy proper- 
ties (including stellar mass, mass-to-light ratio, star forma- 
tion rate, redshift). Wilkins et al. in prep uses a Markov 
Chain Monte-Carlo likelihood exploration algorithm with 
an adaptive Metropolis proposal using the PyMC pack- 
aga*j (BayesME, Zuntz & Wilkins in prep) assuming, for 
the most part, a library built using simple exponentially 
decaying star formation histories (ip ex e~ t ' T ) and a fixed 
metallicity. In the default scenario of Wilkins et al. in prep 
(synthetic galaxies at z = 0.5 with the redshift known a pri- 
ori, broadband observed frame UV-optical-NIR photometry 
with uniform noise where S/N(-H/i60u>)= 100, the l/-band 
mass-to-light ratio was recovered with an accuracy of 0.05 
(where the accuracy is defined as the median of the resid- 
ual distribution R, where R — (log 10 Tv) ~ log 10 F^ uo ), and 
uncertainty 0.05 (where the this is defined as half of the 
15.9 th -84.1 th percentile range of the residual distribution). 
Decreasing the number of observed bands, decreasing the 
S/N or increasing the redshift all result in the uncertainty 
growing. 

The optimal colour relation at z = 0.5 for the HST 
filter set is i/775™ — Zfssoip with pi = 2.37, P2 = —0.51 and 
o R = 0.05. A S/N of 100 in both 1/775™ and 2/ssozp suggests 
o\ r = 2 x (1.0875/100) 2 ~ 0.00012 and thus the uncertainty 
on oigT, in the absence of dust is, using Eq. |13| oi g r — 0.06. 
Degrading the S/N to 30 yields oi gT ~ 0.10. 



https://github.com/pymc-devs/pymc 
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Table 1. The optimum observed- frame colour for estimating the stellar mass-to-light ratio for various rcdshifts assuming the GAMA 
(top) and HST (bottom) filter sets. 1 the observed frame colour chosen to measure the mass-to-light ratio. 2 the rest-frame wavelength 
range. 2 Pearson's correlation coefficient between log 10 IV an d the colour (m a — mj,). 3 The standard deviation of the intrinsic residual 
distribution. 4 Par ameters describing the best fit using simple linear regression (see Equation pi. 5 The gradient of the dust vector d 
(defined in Eq. fill. 6 The factor relating the colour excess Eb — v arid the systematic offset to log 10 T (defined in Eq,|l2|. 
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It thus appears that a single-colour diagnostic is capa- 
ble, despite the availability of significantly less information, 
of recovering the stellar mass-to-light ratio with a similar 
uncertainty to this implementation of SED fitting. While 
this may seem counter-intuitive it reflects the fact that in 
this case improving the sophistication of the fitted model 
does more for the quality of the results than simply adding 
more data. 

A simple intuitive example of this would be an imple- 
mentation of SED fitting which assumes all galaxies have a 
constant star formation history. Such a model is clearly not 
a good description of the real Universe given observations 
of passively evolving galaxies. As such the model is unlikely 
to provide an accurate estimate of the mass-to-light ratio 
(of passively evolving galaxies) . In such a case adding more 
data is unlikely to improve the accuracy/precision of the 
fit and a single-colour diagnostic (based on a more sophisti- 
cated model) would unambiguously provide a more accurate 
diagnostic of the stellar mass-to-light ratio. 



colour diagnostic (cr; 9 r = 0.05 — 0.15) are comparable to 
those obtained using a wider set of observations but sim- 
plistic star formation and metal enrichment histories (e.g. 
exponentially decaying star formation histories). 

For future use we provide details of the optimal choice 
of colour (for two common sets of filters), along with the 
associated parameters necessary to accurately determine the 
stellar mass-to-light ratio at z = — > 1.5. 
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7 CONCLUSIONS 

Using the GALFORM galaxy formation model to produce re- 
alistic galaxy star formation and metal enrichment histo- 
ries we have investigated the correlation between the stellar 
mass-to-light ratio and various observed frame colours. 

We find that several observed-frame colours (typically 
within the resi-frame optical) and the U-band stellar mass- 
to-light ratio are often highly correlated (|r| > 0.8, where 
r is Pearson's correlation coefficient) and that the (V- 
band) mass-to-light ratio can be predicted from the in- 
trinsic colours, using a simple linear relation (log 10 Ty — 
pi(m a — mb)+P2), with an uncertainty of ai g r < 0.07. While 
the addition of dust introduces both a systematic offset and 
increases the uncertainty, if the choice of filters is made care- 
fully these effects are typically small. 

The uncertainties predicted to arise using a single- 
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also provide parameters relating the observed frame colour 
to the observed frame mass-to-light ratio, i.e.: 

log 10 rv b = \og 10 [{M/M e )/(L b /L b {M b = 0.0))] (Al) 
= pi x (m a - m b ) +p 2 . (A2) 

For each redshift/colour the redder filter is chosen as the 
reference band for the luminosity (L b ) as there is typically 
less variation in the mass-to-light ratio. The units of lumi- 
nosity are chosen to be AB centric (as opposed to the more 
familiar solar units in the case of rest-frame V-band mass- 
to-light ratio) and are defined as the luminosity at M b = 0.0 
(i.e. L b (M b — 0.0)). With the exception of the parameter di 
(defined in Eq. Ill, which has to be re-defined as 



di = 0.4 x k b /(k a - k b ), 



(A3) 



(i.e. replacing the kv term with k b ) we can derive the same 
set of parameters (r, an, pi, p2, cfe) that were used to de- 
scribe the rest-frame l/-band mass-to-light ratio - single 
colour relation for the observed frame relation. These pa- 
rameters are given in Table [AT] 



This paper has been typeset from a TpX/ P/IpX file prepared 
by the author. 



APPENDIX A: OBSERVED FRAME 
MASS-TO-LIGHT RATIO 

To avoid the need for a fc-correction (from an observed frame 
luminosity to the rest-frame l/-band) when estimating stel- 
lar masses (as opposed to the stellar mass-to-light ratio) we 
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Table Al. Parameters relating the observed-frame colour (m a — mj,) with the observed frame mass-to-light ratio for the same colours, 
redshifts, and filter sets in Table [I] The observed frame reference band is chosen to be redder filter of each colour (i.e. the fe-band). 
1 the observed frame colour chosen to measure the mass-to-light ratio. 2 Pearson's correlation coefficient between logjg T;, a j, and the 
colour (m a — mj). 3 The stan dard deviation of the intrinsic residual distribution. 4 P aram eters describing the best fit using simple linear 
regression (see Equation |A1b. 5 The gradient of the dust vector d (re-defined in Eq . |A3| for the observed frame b- band ). 6 The factor 
relating the colour excess Eg_y and the systematic offset to log 10 T (defined in Eq,|12[l. 
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